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We demonstrate accurate single-qubit control in an ensemble of atomic qubits trapped in an optical 
lattice. The qubits are driven with microwave radiation, and their dynamics tracked by optical 
probe polarimetry. Real-time diagnostics is crucial to minimize systematic errors and optimize the 
performance of single-qubit gates, leading to fidelities of 0.99 for single-qubit n rotations. We show 
that increased robustness to large, deliberately introduced errors can be achieved through the use of 
composite rotations. However, during normal operation the combination of very small intrinsic errors 
and additional decoherence during the longer pulse sequences precludes any significant performance 
gain in our current experiment. 
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I. INTRODUCTION 

Optically trapped neutral atoms have been widely con- 
sidered as a physical platform for the implementation 
of quantum information processing. Qubits encoded in 
the hyperfine ground manifold offer long coherence times 
and are readily manipulated using the tools of laser spec- 
troscopy, whereas the two-qubit interactions required for 
entangling quantum gates are more difficult to implement 
due to the weak coupling between atoms in the electronic 
ground state. For this reason, two-qubit quantum logic is 
generally regarded as the core challenge for the neutral 
atom platform, and most theoretical and experimental 
efforts have focused on developing viable gate schemes 
based on e. g. cold collisions [I] [3 0] or dipole-dipole 
interactions between Rydberg states [S] [5] [7] [H] . From 
these first-generation experiments it is clear, however, 
that accurate single-qubit manipulation is already criti- 
cal and not always trivial to achieve in a laboratory set- 
ting. In practice, noisy or inaccurate driving fields and 
trap potentials, along with other perturbations, cause 
unpredictable errors in the effective Rabi frequency and 
detuning for the qubit two-level system, and reduce the 
fidelity of rotations on the Bloch sphere. In this situation 
it is useful to explore both how far one can reasonably 
go through a direct effort to eliminate these sources of 
error, and how one can use composite pulse techniques 
that are robust against errors as is routinely done by the 
NMR community [9 [10 , and more recently for qubits 
encoded in trapped ions [H][l2], electron spins [13] and 
superconducting circuits |14j . Ultimately, very high fi- 
delity single-qubit operations can become an important 
resource to combat errors and dephasing across the qubit 
ensemble [TB], through spin echoes and refocusing pulses 
as used in NMR quantum computation [10] [E]. It is 
likely that trapped atom implementations will need to 
adopt a similar strategy. 



For ensembles of atomic qubits trapped in optical lat- 
tices, precise single-qubit control is compromised by spa- 
tial variations in the lattice potential and the associated 
inhomogeneous broadening of the AC Stark shifted qubit 
transition frequency. Furthermore, some collisional gate 
schemes require that a qubit be encoded in magnetic 
field-sensitive states, such as |0) — \F — 3, my = 3) 
and |1) = \F = 4, m/ = 4) in the case of the Cs atom 
qubits used in our work here, in order to permit state- 
dependent motion in the optical potential and/or to sup- 
press spin changing collisions |2J. In that case imperfect 
control over the magnetic field can easily become an im- 
portant source of error. In this letter we demonstrate 
that microwave-driven single-qubit gates with fidelities as 
high as 0.99 can be achieved in parallel on ~ 10 7 atomic 
qubits, by carefully optimizing the spatial and temporal 
uniformity of the optical lattice and of the applied and 
ambient magnetic fields. In doing so, a critical element 
has been the development of a non-perturbing method to 
observe qubit ensemble dynamics in real time and with 
good signal-to-noise ratio (SNR), based on optical probe 
polarimetry. We further show that composite pulse tech- 
niques can significantly increase robustness against pulse 
imperfections and ensemble inhomogencities. One pulse 
sequence in particular, known as a rotary echo, greatly re- 
duces sensitivity to these errors and allows easy observa- 
tion of the coherence time in our system. We expect that 
increased microwave irradiance, in combination with de- 
creased scattering of lattice photons and the use of more 
advanced pulse techniques [T7], will allow us to improve 
the fidelity of single-qubit operations by at least another 
order of magnitude. 

The remainder of this paper is organized as follows. 
In Sec. II we describe our experimental apparatus, in- 
cluding our programmable microwave source and optical 
probe setup. Sec. Ill discusses Rabi oscillations and the 
fidelity of simple rotations, Sec. IV describes the use of 



rotary echoes to extract coherence times, and Sec. V dis- 
cusses the use of composite pulses to improve gate fideli- 
ties in the presence of errors. We summarize our findings 
and discuss the outlook for robust control of neutral atom 
qubits in Sec. VI. 



II. EXPERIMENTAL SETUP 




Our basic experimental setup is shown in Fig. [TJa). 
We begin by preparing a sample of ~ 10 7 laser cooled Cs 
atoms in a cloud with a diameter of ~ 0.8mm. The atoms 
are loaded into a three-dimensional optical lattice tuned 
130GHz to the blue of the 65 1/2 (F = 4) -> 6P 3/2 (F' = 
5) transition of the Cs D2 resonance line, by ramping up 
the lattice potential while the atoms are being simultane- 
ously cooled by a near-resonant 3D optical molasses, con- 
tinuing to cool the atoms for a few ms while in the lattice, 
and then extinguishing the optical molasses beams. The 
optical lattice is formed by three independent, linearly 
polarized standing waves whose frequencies differ by at 
least 10MHz in order to eliminate interference between 
them. Typical lattice depths are Ul ~ IOOEr , and typ- 
ical vibrational frequencies in the lattice potential wells 
are lo l ~ 20E R /h ~ 2tt x 40kHz, where E R = (hk) 2 /2M 
is the single photon recoil energy. Once trapped in the 
lattice the atoms are optically pumped, producing a 90% 
population of the |0) = \F = 3,m/ = 3) qubit state. No 
additional cooling of the center-of-mass motion is per- 
formed, and measurements of the kinetic temperature 
indicate a mean vibrational excitation n ~ 1 for each 
degree of freedom. At this point the qubits can either 
be kept in the lattice, or released into free fall to elim- 
inate lattice-induced light shifts of the qubit transition 
frequency. 

Single-qubit rotations are implemented in parallel 
across the entire ensemble by illuminating the atoms 
with microwave radiation. We use an inexpensive pro- 
grammable microwave source consisting of an HP8672A 
synthesized signal generator running at 9.2GHz, mixed 
with a ~ 30MHz signal from an SRS DS345 synthe- 
sized arbitrary waveform generator. The DS345 provides 
both frequency fine-tuning and arbitrary phase modu- 
lation (0.001° resolution) under control of an uploaded 
computer generated waveform, and, if desired, can be 
amplitude modulated under control of an arbitrary ana- 
log waveform generated by, e. g., a second DS345. By 
driving the qubits with a modulation sideband we thus 
achieve the precise, programmable phase control required 
to implement composite pulses. A digital switch inserted 
between the 9.2GHz local oscillator and the mixer allows 
the mixer output to be turned on/off with a transition 
time of ~ 3ns. The output from the mixer is ampli- 
fied to ~ 2W by a solid state power amplifier, and radi- 
ated by a gain horn located approximately 20cm from the 
atomic sample. Due to reflections from various metallic 
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FIG. 1: (Color online) (a) Experimental setup, showing the 
optical lattice and the configuration used for optical probe 
polarimetry. (b) Birefringent phase shift ip seen by the probe 
(solid line), and the probe induced differential light shift AU 
of the qubit states (dashed line), as function of the probe 
detuning from the 6S l i/ 2 (F = 4) — » 6P 1 / 2 (F' = 4) transition 
of the Dl line. 



surfaces in our apparatus we have only moderate control 
over the microwave polarization, but transitions to hy- 
pcrfine states outside the qubit basis arc Zccman shifted 
out of resonance by a ~ 400mG bias magnetic field and 
effectively suppressed. With this setup we can achieve 
resonant Rabi frequencies as high as 40kHz when driv- 
ing the \F = 3,mf = 3) — » \F = 4,rrif = 4) transi- 
tion, which corresponds to the largest magnetic dipole 
moment. 

High fidelity rotations of the qubit pseudospin require 
precise knowledge of the resonant Rabi frequency x an d 
the effective detuning A of the microwave driving field 
relative to the qubit transition frequency. In our setup 
X varies slightly across the ensemble, due in part to the 
spatially inhomogeneous radiation pattern from the gain 
horn, and in part due to interference with microwave re- 
flections. Variations in A are caused by inhomogeneous 
broadening of the qubit transition frequency, which is 
light shifted by the spatially inhomogeneous lattice light 
field, and Zceman shifted by the applied and ambient 
magnetic fields. For ensembles in far detuned lattices, 
the inhomogeneous spreads in x an d A easily become 
the most important factors limiting the fidelity of single- 
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qubit rotations, and they must be carefully minimized to 
achieve good gate performance. We accomplish this by 
driving Rabi oscillations between the qubit states, and by 
carefully aligning our lattice beams and aiming our mi- 
crowave horn to maximize both the Rabi frequency and 
the decay time for the Rabi oscillation amplitude. In do- 
ing so, an essential requirement is to have good diagnostic 
tools so that adjustments can be made and improvements 
detected as near as possible in real time. 

We perform real-time, minimally perturbing measure- 
ments on our qubit ensemble by optical probe polarime- 
try, as described in detail in |18j . In this measurement 
scheme we choose the atomic quantization axis along the 
(vertical) z-axis, pass a weak, off-resonance probe beam 
through the ensemble along the y-axis, and choose the 
probe polarization to be linear and oriented at 45° in 
the x-z plane. The probe frequency is tuned in-between 
the 6S 1/2 (F = 4) -> 6P 1/2 (F' = 3,4) transitions of the 
Dl line, where the tensor polarizability is substantial for 
atoms in the |1) = \F = 4, to/ = 4) qubit state. In 
this situation the ensemble becomes birefringent by an 
amount proportional to the population of the logical- 11) 
state, and the probe polarization acquires a small degree 
of ellipticity, which is easily detected with a shot-noise 
limited polarimeter. As demonstrated previously for the 
states involved in the atomic clock transition [18], any 
probe-induced, spatially inhomogeneous light shifts of 
the qubit transition frequency can be minimized by tun- 
ing the probe frequency to the appropriate magic value, 
in this case 135 natural linewidths (633MHz) to the red 
of the F = 4 — > F' = 4 transition (Fig.JTJb)). In principle 
one can also detect the Faraday rotation caused by qubits 
in the logical-|l) state, a variation of probe polarimetry 
which is perhaps more familiar from its widespread use 
in magnetometry |19j and in measurements on hyperfine 
spin ensembles [20] [21] [22] • However, for our qubit en- 
coding the Faraday signal almost vanishes at the magic 
probe frequency and is therefore not as useful. Over- 
all, our measurement scheme provides high bandwidth 
access to the qubit dynamics, with little or no pertur- 
bation except for optical pumping caused by the scat- 
tering of probe photons. Even for our relatively dilute 
atomic samples, the SNR's needed to align the lattice 
and optimize our apparatus for minimal qubit transition 
frequency broadening can be achieved in seconds. Much 
higher SNR's are needed to measure the fidelity of highly 
accurate qubit rotations, but still require no more than a 
few minutes averaging. As a consistency check or when 
more detailed information is needed, we complement the 
optical probe measurement with Stern Gerlach measure- 
ments [33] performed at discrete moments in time. While 
time consuming, these provide access to the individual 
populations of all magnetic sublcvcls in the entire ground 
hyperfine manifold, and provide an absolute calibration 
for the initial preparation of the qubits by optical pump- 
ing. 



III. RABI OSCILLATIONS AND SIMPLE 
ROTATIONS 

An arbitrary qubit rotation R u (9) can be specified by 
the axis u and angle 9 of rotation. The fidelity of a lab- 
oratory implementation of R u (9) can be tested by ini- 
tializing an ensemble of qubits in a state \ipi), applying 
the microwave control pulse, and measuring the absolute 
square of the overlap, 9 = \(ip exp \R u (9)\ipi)\ 2 , between 
the target state Ru(9)\i>i) and the state \ip e x P ) produced 
in the experiment. In general, one should sample the fi- 
delity for a distribution of initial states, but for simple (as 
opposed to composite) rotations the fidelity depends only 
on the errors in u and 9, and we can choose = |0) 
without loss of generality. Specializing further to rota- 
tions by an angle 9 — tt around an axis in the equatorial 
plane of the Bloch sphere (a 7r-gate) , the fidelity can be 
determined simply by measuring the populations of the 
logical- 10) or -|1) states. In practice, accurate determi- 
nation of near-unit fidelities is best done by measuring 
the cumulative error in many successive applications of 
the rotation. Experimentally, this is equivalent to driv- 
ing Rabi oscillations and observing their decay due to 
decoherence and dephasing across the ensemble. 

Our atomic qubits are embedded within a larger 
ground hyperfine manifold, and can be coupled to states 
outside the logical space either by the microwave driving 
field (see below) or by optical pumping due to scatter- 
ing of photons from the probe and optical lattice light 
fields. It is straightforward to solve the master equation 
for a microwave driven atom in the presence of optical 
pumping; we have done so for parameters typical of our 
experiment and have found that the qubit pseudospin 
can be approximated to an excellent degree as a lossy 
2-lcvel system. This allows us to replace the numerically 
intensive master equation model with a much simpler 
effective model based on the Torrey solutions to the fa- 
miliar Bloch equations |24j |25j . In our experiments the 
microwave resonant Rabi frequency is always much larger 
than the inverse of the longitudinal (Ti ) and transverse 
(T 2 ) homogeneous lifetimes, so that \ ^> 1/^2 — l/^i (the 
strong driving limit of the Torrey solutions) . The detun- 
ing of the microwave field relative to the qubit transition 
frequency is also relatively small, A ~ 10 _1 . With the 
initial condition \ipi) — |0), the Torrey solution for the 
expectation value of the 3-component of the qubit pseu- 
dospin (the inversion) is then 

w ® = -^ ex p^ f71t +iJ cos ( m ) ex p~ 71t + 2 ^ (!) 

where fl = \/ x 2 + A 2 and 71 = 3/(4Ti), and where we 
have set T' 2 = 2Ti as appropriate for a closed two-level 
system subject only to optical pumping. If the total pop- 
ulation of the two- level system decays at a rate 72, then 
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the logical- 1 1) population oscillates in time as 

ni(t) = ^>W + i]e" 72t 



1A= 

2 n 2 



-(§71+72)* 



"72* 



2 ft 2 



cos(fii)e- (7l+72)t , 



(2) 
(3) 



where the last approximation introduces errors of order 
10~ 3 for the parameters typical of our experiments. 

Inhomogeneous effects occur due to variations in the 
microwave resonant Rabi frequency and/or detuning. We 
therefore define an ensemble averaged population 



/oo />oo 
d x / dAP(x,A)n!(t), 
-oo J — oo 



(4) 



where the distribution P(x, A) is Gaussian, with mean 
values Xoi A and standard deviations Sx, SA for the 
resonant Rabi frequency and detuning, respectively. In 
the end, our expected polarimetry signal is of the form 



s(t) = sMt) + s 



(5) 



where the signal amplitude Si and offset So depend on 
probe power and balancing of the polarimeter. The func- 
tion S(t) is used to fit our experimental data sets, with 
free parameters XOi Ao, Sx, 8A, Si, and So. The homo- 
geneous decay rates 71 and 72 are not included as free 
parameters in the fit, but instead determined indepen- 
dently in a separate spin-echo experiment as described 
in Sec. IV. 

Figure [2] shows typical examples of our polarimetry 
signal for atoms in free fall, and for atoms trapped in 
the optical lattice. Two features arc immediately appar- 
ent in both data sets: the oscillation amplitude decays 
in a few milliseconds, corresponding to several tens of 
Rabi periods, and the mean of the oscillation decays on 
a much slower timescale. This is characteristic of situa- 
tions where dephasing plays a significant role. The figure 
also shows fits of the form S(t) to each signal. In the case 
of free atoms such fits imply a microwave irradiance in- 
homogeneity across the atomic sample of Sx/xo = 0.3%, 
and a detuning inhomogeneity, SA/xo — 3.3%, presum- 
ably caused by spatial variations in our bias magnetic 
field. For atoms trapped in the optical lattice we find 
$x/Xo = 0.3% and SA/xo — 7.3%, with the increase in 
(5A reflecting the additional inhomogeneous broadening 
from spatial variations in the lattice light shift. These 
values for free and trapped atoms represent averages over 
a number of measurements taken under nominally iden- 
tical conditions but at widely different times. 

Based on the parameters xoj Ao, 5x, and 5A extracted 
from high-quality fits to Rabi oscillation data, and on the 
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FIG. 2: (Color online) Polarimetry signal S(t) measured dur- 
ing qubit Rabi oscillation, for atoms in free fall (top) and 
atoms trapped in the optical lattice (bottom). The dark line 
is experimental data, and the underlying light line is the result 
of a fit of the form given in Eq. |Hj. Inserts show a close-up 
of data points (filled circles) and fit for a representative time 
interval. 



homogeneous decay rates 71 and 72 determined from spin 
echoes, we can accurately model the ensemble averaged 
population fl 1 (t) of the logical-|l) state. This in turn 
allows us to obtain a good estimate of the ensemble av- 
eraged fidelity j^. for a perfectly timed 7r-gate, defined 
as a square microwave pulse of duration r = tt/xo a P~ 
plied across the entire ensemble. Note that in doing so, 
the relevant decay rates are those obtained in the pres- 
ence of the optical lattice only, since the probe beam 
would not be present in experiments aimed at quantum 
information processing. We have, in this fashion, found a 
7r-gate fidelity = 0.990(3) for atoms in the optical lat- 
tice, averaged over 16 independent measurements taken 
over a 3 month period. This value is a measure of the 
limitations imposed by working with qubits in spatially 
inhomogeneous lattice and microwave fields. 

In practice, the fidelity of a gate operation is further 
reduced if we cannot preset the mean Rabi frequency xo 
and detuning Ao with sufficient accuracy. From our data 
set we find statistical fluctuations of 1% rms in the fit- 
ted value for the mean Rabi frequency xo, relative to our 
target value of 27.78kHz which corresponds to an 18.0/is 
7r-gate time. In the same data set the rms fluctuations in 
the mean detuning were Aq/xo ~ 4.8%. We can use the 
fit parameters for each individual measurement to calcu- 
late IIi(18.0/is), which is the ensemble-averaged fidelity 
at the time intended to produce a 7r-gate. Averaged over 
our 16 independent measurements, this reduces the fi- 
delity by a statistically insignificant 1 x 10 -4 . 

In our discussion so far, we have ignored the possibility 
of microwave coupling to states outside the qubit logical 
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space. For our encoding such coupling occurs only be- 
tween the |F = 3, itif = 3) logical- 10) state and the states 
\F = 4, nip = 2, 3) . These transitions arc Zccman shifted 
out of resonance by the applied bias magnetic field, and 
further suppressed by maximizing the cr + -polarized com- 
ponent of the microwave field. For a quantitative es- 
timate, it is straightforward to tune our microwave fre- 
quency to be resonant with these transitions and to deter- 
mine the resonant Rabi frequencies in each case. For the 
data sets discussed here, we thus found Rabi frequencies 
of 3kHz for both transitions, and detunings of 130kHz 
and 260kHz respectively. During a 7r-gate on the qubit 
this translates into a combined population of 7 x in 
the two states, and thus a reduction in the 7r-gate fidelity 
of less than 1 x 1CP 3 . This is well below the statistical 
uncertainty of our present experiment, but if necessary 
the population leakage can be further suppressed with a 
stronger bias magnetic field. 

IV. ROTARY ECHOES AND COHERENCE 
TIMES 

As the preceding discussion makes clear, the decay of 
qubit Rabi oscillations in our optical lattice is dominated 
by dephasing due to inhomogeneities in the microwave 
irradiance and light shift across the ensemble. In this 
situation the underlying coherence times cannot be re- 
liably determined by fitting experimental measurements 
of Rabi oscillations with the model given in To cir- 
cumvent this problem we use a technique originally de- 
veloped by the NMR community, rotary echoes [TU], to 
suppress dephasing and obtain independent estimates for 
the rates 71 and 72. The simplest rotary echo sequence 
consists of repeated application of the composite rotation 
R u (8)R u i(9), where the axes of rotation are controlled 
through the phase tp of the driving field. For given Rabi 
frequency and detuning we have 

u = ^[cos(y>)l + sh%)2] + ^3 (6) 

In the near resonance case, A <C 0, a phase change 
if — > ip + 7r will nearly invert the axis of rotation, u' « u, 
and the composite rotation remains close to the identity. 
Thus, if the qubit is initialized in logical- 10) and is an 
integer multiple of 7r, the rotary echo sequence will pro- 
duce a polarimetry signal that is indistinguishable from 
standard Rabi oscillations with greatly suppressed inho- 
mogeneity in \- hi general, rotary echoes dephase as 
rapidly due to inhomogeneity in A as do Rabi oscilla- 
tions. Remarkably, however, the special case 8 — 2tt is 
robust against variations in \ and A, and a rotary-27r 
echo sequence looks like Rabi oscillations with greatly 
suppressed inhomogeneity in both. This is easy to verify 
through numerical simulation of rotary-27r echoes, which 
for inhomogeneities typical of our experiment show no 
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FIG. 3: Polarimetry signal S(t) during qubit Rabi oscillation 
(top) and a rotary-27r spin-echo sequence. An artificially large 
inhomogeneity in the microwave detuning A has been added 
to demonstrate the robustness of the rotary-27r echo. 

sign of dephasing on timescales far beyond the longest 
coherence times that we have observed in the laboratory. 
Figure [3] is an experimental illustration of this robust- 
ness, comparing polarimetry signals from Rabi oscilla- 
tions and a rotary-27r echo in the presence of an artifi- 
cially increased inhomogeneity in A. 

The analysis of qubit Rabi oscillations in Sec. Ill re- 
quires that we determine the rates 71 and 72 in the pres- 
ence of a probe beam of identical intensity and frequency. 
This is easily achieved by recording the measurement 
from Rabi oscillation and a rotary-27r echo in quick suc- 
cession, and fitting the latter by setting S\ — SA = in 
our model for Rabi oscillation. To estimate gate fideli- 
ties we also need to determine these decay rates in the 
limit where photon scattering from the probe becomes 
negligible. Figure |3Ja) shows rotary-27r echo signals for a 
few probe powers, clearly illustrating the decrease in sig- 
nal strength and increase in coherence time as the probe 
scattering rate is reduced. Figure |4^b) shows the time 
constant Td for the decay of the oscillating part of the 
echo signal, as a function of the scattering time (the in- 
verse probe scattering rate) t s = l/j s . As one would 
expect, Td is roughly proportional to r s for short scatter- 
ing times, and saturates for long scattering times. Also 
shown is a fit of the form Td = 1/ (a7s + b), which for this 
data set gives an asymptotic value Td ~ 10.6ms that pre- 
sumably reflects the coherence time for qubits in the op- 
tical lattice. Note that this is only a representative data 
set; over time and for nominally identical conditions we 
have observed asymptotes in the range 5.5ms— 15ms. For 
the data analysis in Sec. Ill we therefore used the most 
conservative estimate, Td = 5.5ms. For technical reasons 
(low-frequency noise and weak polarimetry signals at low 
7 S ) it is more challenging to get an accurate estimate for 
the decay of the constant part of the echo signal. We 
have generally found this decay time to be ~ 2 x Td, and 
therefore 71 « 72 ~ l/(2rd) according to Eq. ([3]), and 
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FIG. 4: (a) Polarimetry signals S(t) measured during rotary- 
2tt spin-echo sequences. The probe power is reduced by a 
factor of ~ 40 from top to bottom, to demonstrate the longer 
coherence times and lower signal levels available at lower rates 
of probe photon scattering, (b) Coherence time Td versus 
probe photon scattering time r s . 



we have assumed this throughout our data analysis. This 
observation is also consistent with a full master equation 
model of the dynamics of the entire &Si/ 2 ground mani- 
fold in the presence of optical pumping and a microwave 
drive. 

It is instructive to consider what physical mechanisms 
might limit the qubit coherence time. First, for our en- 
coding the qubit transition frequency is sensitive to tem- 
poral fluctuations in the lattice depth and in the back- 
ground magnetic field. As long as the power spectrum 
of these fluctuations lies well below the Rabi frequency, 
driving continuous Rabi oscillations or rotary echoes is 
equivalent to applying a sequence of refocusing pulses, 
and one expects to see longer coherence times than in, e. 
g., a Ramsey interrogation [IS]. Our experiment does in- 
deed show such an increase. Second, qubits decohere due 
to optical pumping caused by the scattering of photons 
from the optical lattice. For our relatively deep and near- 
resonance optical lattice we calculate a scattering time of 
~ 100ms, which is too large to account for the ~ 10ms 
coherence times seen in our experiment. Further evidence 



of the minimal role played by lattice photon scattering 
comes from the fact that the coherence time remains ef- 
fectively constant when the scattering rate is changed by 
as much as a factor of two. Finally, we note that phase 
noise in the microwave source will cause the axis of rota- 
tion to fluctuate and thus degrade the rotary echo. If the 
phase noise spectrum is flat we expect the echo signal to 
degrade after a roughly constant number of echo pulses, 
resulting in longer coherence times when the microwave 
Rabi frequency is reduced. The qualitative behavior of 
our experiment suggests that this may be an important 
contributing factor in limiting the coherence time, but 
our attempts to directly measure the phase noise and 
model its influence on the qubit dynamics have so far 
been unsuccessful. 



V. COMPOSITE PULSES 

Composite rotations are a powerful tool to overcome 
the effect of inhomogeneous errors, not just in spin-echo 
and refocusing sequences, but also in the implementation 
of quantum gates. Cummins [5] describes three families 
of universal pulse sequences which are designed to con- 
fer robustness in different scenarios. The term univer- 
sal here refers to a pulse sequence that can implement 
any desired unitary transformation, i. e. any rotation 
about a given axis by a given angle, with a similar de- 
gree of error compensation for all initial states. This 
contrasts with pulse sequences designed to be robust only 
for the transfer of a specific input state to a specific out- 
put state, e. g. the familiar 3-pulse refocusing sequence 

i^/2)#2(V>-Ri(>/2)- 

For our system the dominant factor limiting gate fi- 
delity is the variation of effective microwave detuning 
across the qubit ensemble. The so-called CORPSE 
(Compensation for Off- Resonance with a Pulse Se- 
quence) family is a three-pulse sequence that confers ro- 
bustness against these detuning errors. If, for example, 
our goal is to perform a perfect it rotation around the 1- 
axis, then we can approximate this through the CORPSE 
sequence 



i?u(^/3)i?u'(57r/3)i? u (77r/3) « Ri(tt), 



(7) 



where u = (|x|/fi)l+(A/n)3, u' = -(| X |/f2)l+(A/0)3. 
Figure [5] shows the corresponding trajectory on the Bloch 
sphere, as well as the evolution of the logical- 11) pop- 
ulation during the pulse. To leading order in the de- 
tuning, / = A/xo *C 1, the fidelity of a single square 
pulse is J*" « 1 - f 2 /2. By comparison, a CORPSE 6- 
pulsc achieves & ss 1 — a/ 4 , where a = 6.5 x 10~ 3 for 
O = 7r, at the cost of increased sensitivity to errors in 
the rotation angle. Similarly, the SCROFULOUS (Short 
Composite Rotation For Undoing Length Over and Un- 
der Shoot) family is a three-pulse sequence that confers 
robustness against errors in the rotation angle at the cost 
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FIG. 5: (Color online) (a) Trajectories on the Bloch sphere 
during the three parts of a composite CORPSE 7r-pulse, in 
the presence of a detuning A = 0.3x- (b) Evolution of the 
logical- 1 1} population during the pulse. 



of increased sensitivity to detuning errors, and the BB1 
(BroadBand Number 1) family is a five-pulse sequence 
that compensates angle errors to sixth order at little or 
no cost in sensitivity to detuning errors. We refer the 
reader to the literature [5] [TO] for further discussion of 
these pulse families, including how to design composite 
rotations by angles other than tt. We note also that new 
numerical techniques enable the brute-force optimization 
of much more complicated time dependent phase varia- 
tions, which can improve the robustness to specific com- 
binations of amplitude and detuning errors, or achieve 
better robustness in the presence of decoherence [17] . 

It is straightforward to check the performance of com- 
posite pulses in our system by preparing our atoms in 
logical- 1 0), apply the appropriate pulse sequence, and 
measure the logical- 11) population. Figure [6] shows the fi- 
delity of plain pulse, CORPSE, SCROFULOUS and BB1 
7r-gates for large, intentionally applied errors. Each ex- 
perimental data set closely follows the predictions of our 
theoretical model, and illustrates how a given pulse fam- 
ily significantly extends the range of detuning or angle 
error over which one can achieve good fidelity. It is worth 
noting, however, that for SCROFULOUS and BB1 pulses 
the peak fidelity falls somewhat below that of a plain 
pulse. This reduction is due to a combination of fac- 
tors: these composite pulses protect against angle errors, 
which are insignificant when the intentional error is near 
zero, the SCROFULOUS pulse increases sensitivity to 
the intrinsic detuning spread in our experiment, and both 



take longer to execute (for a given Rabi frequency) and 
therefore subject qubits to additional decoherence during 
the gate operation. In our case, the gain from increased 
robustness is more than offset by the loss from decoher- 
ence, and the overall fidelity goes down. The tradeoff is 
more favorable for CORPSE pulses, where there is some 
gain from increased robustness to the intrinsic detuning 
errors that dominate in our experiment. 

High fidelity CORPSE 7r-pulses can be evaluated in 
much the same way as plain tt- pulses: by successive ap- 
plication and measurement of the logical- 10) and -|1) pop- 
ulations. Figure [7j a) compares polarimetry signals for a 
rotary-2-7r echo, and for successive CORPSE pulses. Note 
the similar decay times for the two pulse sequences, as 
expected when detuning inhomogeneity dominates. In 
principle the CORPSE data can be analyzed using an 
approach analogous to Sec. Ill, but in practice it is diffi- 
cult to fit this rather complicated polarimetry signal. In- 
stead, we read off the logical- 10) and -|1) populations at 
the beginning and end of each CORPSE pulse, and com- 
pare these to the predictions of a full model based on the 
optical Bloch equations, using the inhomogeneities and 
relaxation rates determined from Rabi oscillations and 
rotary-27r echoes performed under identical conditions. 
Figure [7jb) shows a typical data set for which these val- 
ues agree quite closely. This provides confidence that our 
model correctly describes the performance of CORPSE 
pulses in our experiment, and allows us to estimate a 
CORPSE fidelity ^corpse = 0.992. Since the inhomo- 
geneities and relaxation rates for this data set are typical 
of our larger data set, we conclude that CORPSE pulses 
will yield only marginal improvement for our conditions. 
However, in situations where inhomogeneities are larger 
or coherence times longer, it appears likely that CORPSE 
or other composite pulse sequences can significantly im- 
prove the fidelities of single-qubit rotations. 

VI. SUMMARY AND OUTLOOK 

We have encoded qubits in the ground state manifold 
of neutral Cs atoms trapped in an optical lattice, and im- 
plemented accurate single-qubit 7r-gates simultaneously 
across the entire ensemble by driving the atoms with mi- 
crowaves from a phase agile, programmable source. The 
fidelity of these gates was found to depend primarily on 
spatial variations in the qubit transition frequency which 
is Zeeman shifted by applied and background magnetic 
fields, and light-shifted by the presence of the optical lat- 
tice. Spatial variations in the microwave Rabi frequency 
and decoherence from scattering of lattice light played a 
lesser role in limiting the fidelity of simple rotations. Ac- 
curate estimates of the 7r-gate fidelity were obtained from 
quantitative modeling of Rabi oscillation data. When av- 
eraged over many data sets taken over a 3 month period, 
this yielded a value J^- = 0.990(3). Independent esti- 
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FIG. 6: (Color online) Fidelity of a plain and three types of composite 7r-pulses, as function of deliberately applied errors in either 
the average detuning Ao or average resonant Rabi frequency xo- The pulse sequences are (a) CORPSE, (b) SCROFULOUS, and 
(c) BB1. Filled and open circles are experimentally measured fidelities for the composite pulses and plain pulses, respectively. 
Solid lines are theoretical fidelities found by solving the optical Bloch equations, using independently determined inhomogeneities 
and coherence times. 



mates of coherence times were obtained from rotary-27r 
spin echo sequences that effectively compensates for in- 
homogeneities. Finally, we implemented three types of 
composite pulses and tested their ability to compensate 
for large, deliberately applied errors in Rabi frequency or 
detuning. The CORPSE sequence provides increased ro- 
bustness against detuning errors, and is potentially the 
most useful in an atom/lattice system where the plain 
pulse fidelity is limited by detuning inhomogeneity across 
the ensemble. In our experiment the inhomogeneities 
and coherence times are such that the advantage from 
increased robustness is outweighed by additional deco- 
herence during the longer CORPSE sequence, and we 
were not able to achieve any significant gain in fidelity 
However, this situation may not be representative, and 
optical lattice experiments elsewhere could offer tradeoffs 
that favor the use of composite pulse techniques. 

The single-qubit gate fidelities demonstrated here are 
probably good enough that they will not limit neutral 
atom quantum information processing in the near fu- 
ture. Nevertheless, single-qubit control being the sim- 
plest task on this particular platform, it is important to 
consider how one might improve the performance signif- 
icantly. The most straightforward approach is probably 
to increase the microwave Rabi frequency, by using a 
more powerful microwave source and/or moving the gain 
horn closer to the atomic sample. Increasing the Rabi 
frequency ten-fold while keeping other parameters un- 
changed should reduce the single-pulse 7r-gate time to 
~ 5.7/is, and improve fidelity by nearly an order of mag- 
nitude, to & = 0.9987. Next, one might attempt to 
decrease inhomogeneities in Rabi frequency and detun- 
ing, most likely by working with a more compact sample 



of atoms, or to increase the coherence time enough for 
the robustness versus decoherence tradeoff to favor the 
use of composite pulses. Much longer coherence times 
should be readily achievable if one chooses a magnetic 
field insensitive qubit encoding, employs far detuned op- 
tical traps with lower scattering rates, and improves the 
phase noise of the microwave source. Furthermore, the 
use of relaxation optimized pulse sequences has the po- 
tential to significantly reduce the decoherence incurred 
for a desired level of robustness [T7] . 

A key challenge for quantum information processing in 
optical lattices is to develop methods to address and ma- 
nipulate individual atomic qubits. In large period opti- 
cal lattices such addressable gates might be implemented 
with optical Raman transitions and focused laser beams 
|26j . in a manner similar to that used in ion traps |27j . 
Alternatively, one might use magnetic field gradients [S5J , 
or the light shift in a tightly focused laser beam [33], to 
map qubit position into a shift in transition frequency, 
and use microwaves to selectively manipulate an individ- 
ual qubit or a group of qubits. In practice it will be 
challenging to accurately position a focused laser beam 
relative to an optical lattice, leading to significant uncer- 
tainty in the frequency shift associated with a given qubit 
position. In that case it becomes essential to use ad- 
vanced composite pulse techniques that ensure high gate 
fidelity for the targeted qubit, and simultaneous suppres- 
sion of spurious rotations of the neighboring qubits. The 
CORPSE response of Fig. |6ja) serves as a rough illustra- 
tion of the desired behavior: it offers high fidelity for a 
limited range of detunings, and steep roll-off outside the 
acceptance bandwidth. However, the CORPSE sequence 
is not optimized to suppress qubit response everywhere 
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FIG. 7: (a) Polarimetry signals S(t) measured during 
a rotary-27r spin-echo sequence (top) and a sequence of 
CORPSE 7r-pulses. The nearly identical decay times demon- 
strate that the CORSE pulse successfully compensates for de- 
tuning errors, with little or no cost incurred from the increased 
sensitivity to angle errors, (b) Population of the logical-jl) 
state after an integral number of CORPSE 7r-pulses. Solid 
circles are experimental data extracted from the signal in (a). 
Open circles are theoretical values found by solving the op- 
tical Bloch equations, using independently determined inho- 
mogeneities and coherence times. 



outside the acceptance bandwidth, and is therefore un- 
suitable for robust addressing. Preliminary work suggests 
that numerical optimization of more complex pulse phase 
modulation can be used to design composite pulses with 
nearly ideal response [3D] ■ We hope to explore the use of 
numerically optimized composite pulses in future exper- 
iments aimed at higher fidelity and addressable singlc- 
qubit gates. 
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